Impact of a viscous liquid drop 
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We simulate the impact of a viscous liquid drop onto a smooth dry solid surface. As in experiments, 
when ambient air effects are negligible, impact flattens the falling drop without producing a splash. 
The no-slip boundary condition at the wall produces a boundary layer inside the liquid. Later, the 
flattening surface of the drop traces out the boundary layer. As a result, the eventual shape of 
the drop is a "pancake" of uniform thickness except at the rim, where surface tension effects are 
significant. The thickness of the pancake is simply the height where the drop surface first collides 
with the boundary layer. 



The impact of a liquid drop onto a dry solid sur- 
face lies at the heart of many important technologi- 



cal processes 
spray iii 
bustion chamber 




from the application of a thermal 
to atomization of fuel in a com- 
lH. Recent experiments re- 
vealed the splash formed when a low- viscosity liquid, such 
as water or ethanol, first collides with a dry smooth wall 
at several m/s owes its existence entirely to the presence 
of air 1^, 14, 15, T^]- These results are motivating new 
studies on the large-scale deformations created by im- 
pact when air effects are absent as well as how a splash 
forms [13, H E, 0, [m, m, [13 . 

Here we focus on the impact of a viscous liquid drop 
when air effects are absent. Recent experiments show 
that reducing the ambient gas pressure also suppresses 
the splash of a silicone oil drop. However, the form of 
the splash is very different. While the splash from a low- 
viscosity liquid develops within a few 10 /is of impact, 
the splash from a silicone oil develops slowly, becoming 
evident only after most of the liquid drop has fallen and 
flattened into a thin pancake 0, [25[. We use an axi- 
symmetric Volume-of-Fluid (VOF) code to simulate the 
impact at reduced ambient pressure [2^, H^. Our 
results show that a boundary layer, corresponding to a 
thin region where the radial flow created by impact ad- 
justs to the no-slip condition at the wall, is created by 
the impact. The boundary layer has uniform thickness. 
As impact nears its end, the drop surface flattens onto 
the boundary layer, evolving into a pancake of uniform 
thickness. 

The Volume-of-Fluid simulation solves the Navier- 
Stokes equations, together with constraint of incompress- 
ibility, for both the liquid interior and the gas exterior at 
reduced pressure. Physically appropriate boundary con- 
ditions, in particular the Laplace pressure jump across 
the surface due to surface tension, are enforced [26]. In a 
typical run, an initially spherical liquid drop with radius 
a collides with a dry, smooth solid surface at an impact 
speed Uq oi several m/s. The ambient air density pg 
is kept so small that 2-fold changes in in the value of pg 
have little effect on the liquid dynamics. The bottom sur- 



face of the liquid drop is not broken upon impact, which 
corresponds to maintaining an apparent contact angle of 
180° where the liquid rim meets the wall (Fig. [2^) [2^. 
We require that the velocity field inside the drop satisfies 
both the no-flux and no-slip boundary conditions at the 
solid wall. 

Figure 1 plots successive drop surface shapes (out- 
lined in white) against snapshots from an experimen t by 
Driscoll & Nagel under the same impact conditions [2^. 
We have rescaled time by the impact time-scale r = 
a/Uo. The simulation agrees very well with the exper- 
iment. After the drop hits the wall, the liquid liquid in- 
side the drop is diverted outwards, forming a thin sheet 
which expands radially along the wall (Fig. lb). Later, 
the falling drop flattens, evolving towards a shape resem- 
bling a pancake with a thickened outer rim (Fig. Ic). By 
t = 7.4r, a considerable time after the drop has stopped 
falling, the expanding drop attains its maximum extent 
(Fig. Id). After that point in time, surface tension causes 
the drop in the simulation to retract inwards and reform 
into a spherical shape, a dynamics we do not analyze. 

For this impact, the liquid drop has radius a = 0.16 
cm and the impact speed Uq is 4 m/s. The liquid is a 
low-molecular-weight silicone oil with dynamic viscosity 
PL = 9-4 cP, density pL = 0.94 g/cm^ and surface tension 
cr = 21 dynes/cm. The exterior fluid corresponds to air at 
34 kPa, with density pg = 4.4 x 10~^ g/cm^ and dynamic 
viscosity pg = 1.8x 10~^ cP. The simulated impact occurs 
in a cylindrical tube with radius R = 6a and height H = 
6a. The experiment uses a larger tube. Changing R does 
not change the impact results. In describing the results, 
we use a cylindrical coordinate system, where the tube is 
centered at r = and the wall lies at z = 0. 

To correlate different impacts, we non-dimensionalize 
all the length-scales by a, the velocities hy Uq and the 
time-scales by r. Since neither the tube dimensions nor 
the air properties affect the liquid dynamics reported 
here, the outcomes depend on only two dimensionless 
parameters: the Reynolds number Re = 2pLUoa/ pL and 
the Weber number We = 2pLU'^a/(T. The impact in 
Fig. [1] correspond to Re ^ 1280 and We ^ 2280. 



FIG. 1: Impact of a viscous silicone oil drop at 4 m/s onto a smooth, dry substrate at reduced ambient pressure (34 kPa). 
Surface profiles (white) from simulated impact are overlaid against snap shots from an experiment. From left to right, the 
successive times are t = 0, r, 2r and 7.4r where r = a/[/o is the impact time-scale. Here Re ~ 1280 and We ~ 2280. Photos 
courtesy of Driscoll & Nagel. 



How the thin liquid sheet ejected by impact evolves 
over time is more difficult to characterize in the exper- 
iments. Previous studies 0, ^ proposed a simple es- 
timate for the eventual thickness h of the pancake. The 
idea is that if the only mechanism that can arrest the 
outward radial expansion is viscous dissipation. We then 
assume that the impact energy is dissipated by a radial 
flow of strength [/q in a liquid pancake of thickness h and 
maximal extent i?max- Balancing the dissipation against 
the initial kinetic energy yields 



n 



(1) 



This energy balance, together with volume conservation 
47ra^/3 ^ ^Rmax^^ predicts that the dimensionless pan- 
cake thickness h/a should scale as Re~^^^. This feature 
has been confirmed by previous experiments [31]. Here 
we find that the eventual thickness of the pancake formed 
by impact from our simulation for different values of /7o, 
(J and jiiL are entirely consistent with this estimate (see 
supplementary materials). 

While successful, the scaling estimate for h does not 
tell us when the characteristic pancake thickness h first 
emerges. Nor does it explain why the top surface of 
the drop evolves into a pancake of uniform thickness. 
To address these questions, we examine results from the 
simulation. This time we switch to a reference frame 
X = r — Rh{t). The radial location Rh corresponds to 
a point where the rounded rim joins onto the rest of 
the drop. Within the 0(r) time window, the simula- 
tion shows that Rh{t) ^ y/AaU^t^ so that the reference 
frame x decelerates over time. Figure [2^ displays the 
surface profiles at the leading edge at different moments. 
Initially (t = 0.3r) a thin collar is ejected from a nearly 
spherical drop. As time goes on, a local minimum de- 
velops in the interface profile, separating the drop profile 
into a rounded rim region where surface tension effects 
are important and a downward sloping profile which be- 
comes more and more gently sloped over time. As time 
goes on the profile flattens to the left of the minimum. 
By t = 4.9r, a pancake of uniform thickness is apparent. 
At the rim, surface tension acts to slow the expanding liq- 
uid sheet, causing liquid to accumulate, consistent with 
results from previous studies H, 3^, 3^ . 

The flattening dynamics, however, has not been exam- 



ined previously. To quantify its progress, we define an 
onset time tonset- This is the time when the height of 
the interface at the local minimum first equals the even- 
tual pancake thickness h. In Fig. [2)3 we plot tonset as a 
function of the impact parameters. Non-dimensionalizing 
^onset by the impact time-scale r = a/Uo produces essen- 
tially flat curves with respect to both Re and We. In 
other words, tonset is simply controlled by the kinematics 
of impact, with no apparent dependence on either the 
liquid viscosity or surface tension. As a comparison we 
plot tmax, the time when the drop attains its maximum 
extent and then retracts due to surface tension. This 
quantity has a strong dependence on surface tension, or 
We. Since tonset is considerably shorter than tmax, there 
is still appreciable liquid motion inside the drop when 
the first flattening begins. This suggests that the emer- 
gence of h is not related to whether the kinetic energy 
has been sufficiently dissipated, but instead depends on 
the kinematics. 

To get some insight into the kinematic origin of /i, we 
examine the vorticity field. Prior to impact, the drop is 
falling with a spatially-uniform downward velocity, so the 
vorticity is everywhere inside the drop. After impact, 
the no-flux condition at the wall causes the liquid previ- 
ously falling downwards to be diverted into a radially ex- 
panding flow. This expansion flow speeds up as it moves 
away from the center line, reaches a peak at Rh, and then 
slows down as it enters the rim. This radial expansion 
also adjusts, via viscous effects, to the no-slip boundary 
condition at the solid wall. As a result of this adjust- 
ment, vorticity is generated in the liquid layer nearest to 
the solid wall. At any moment, the amount of vorticity 
generated to ensure zero slip at the wall is dictated by the 
strength of the radial expansion flow. Because the sim- 
ulated impact is axisymmetric, only the azimuthal com- 
ponent of the vorticity is nonzero, i.e. = uj{r^z^t)e0 
where uj = duz/dr — dur/dz. For high Reynolds number 
flows this adjustment takes place inside a narrow bound- 
ary layer. In Fig. [3^ we plot the wall value of the vorticity 
u;o{r^t) [35]. Since the radial outflow is largest near the 
outer edge, the vorticity also has a maximum near Rh. 
As the impact proceeds, the downward fall of the liquid 
drop slows, slowing the expansion and thus reducing the 
magnitude of loq. 

We next outline the spatial extent of the boundary 



3 



layer. This task is complicated by the fact that the ab- 
solute size of the vorticity is strongly correlated with the 
strength of the radial expansion flow. Since the radial 
expansion flow, generated due to the no-fiux condition, 
varies with r, vorticity is generated at different rates at 
different spatial locations. Moreover, at a given loca- 
tion, the vorticity production rate slows over time be- 
cause the radial expansion slows. Thus contours of ab- 
solute vorticity do not provide clear indications for the 
spatial extent of the boundary layer. We side-step this 
complication by normalizing a;(r, the vorticity dis- 
tribution in the bulk of the liquid, by the "wall" value 
iJo{r,t). This essentially strips away variations in the 
vorticity distribution due to the varying speed of the ra- 
dial expansion flow. For simple high Reynolds number 
flows, such as a uniform flow past a solid wall, as well as 
the boundary layer created by a straining flow towards 
a solid wall [36|, this procedure correctly reproduces the 
boundary layer structure that emerges from asymptotic 
analysis. In Fig. 4b we plot contours of the normal- 
ized vorticity distribution. In each snapshot, the solid 
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FIG. 2: Time evolution of ejected liquid sheet, (a) Leading- 
edge evolution in the co-moving frame x = r — Rh{t) . Here 
we reduced the impact speed to [/o = 2 m/s and increased the 
liquid viscosity /hl = 0.2 poise to generate a thicker pancake. 
Other parameters unchanged, (b) Onset time tonset (open 
symbols) and tmax, the time of maximum extent (closed sym- 
bols) as a function of the Weber number (We = plUqO? /a). 
The different symbols correspond to Uq = 2-8 m/s (o), 
cr = 5.25-84 dynes/cm (v), and /i = 10-50 cP (A). Inset 
plots onset time tonset vs. Re = 2aUQ/h'L- 
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FIG. 3: Vorticity evolution after impact, (a) Value of the 
vorticity at the wall as a function of radial distance r [351 . 
From top to bottom, the profiles are taken at t = 0.3t, 0.9t, 
1.6r and 3.0t. (b) The 20%, 50%, and 90% contours of the 
vorticity distribution. 



lines within the liquid drop correspond to contours where 
uj{r,z,t)/uJo{r,t) = 90% (lowest curve), 50% and 20% 
(highest curve). At early times {t = 0.3r), the bound- 
ary layer delineated by the contours is a pancake-shaped 
region inside the liquid drop. Except at the outermost 
edge, the top surface of the liquid drop is widely sepa- 
rated from the boundary layer. As impact proceeds, the 
boundary layer extends radially and thickens slightly, but 
retains its pancake shape. At t = tonset, the top surface 
of the drop collides with the boundary layer. After the 
collision, the surface at the collision location ceases to 
decrease in height. At the same time, the rest of the 
drop surface continues to fall downwards, bringing more 
and more portions of the surface into collision with the 
boundary layer. The result is a "front" that flattens in- 
wards radially, tracing out the pancake shaped boundary 
layer. 

The idea that the eventual pancake thickness h first 
emerges when the top surface of the drop collides with 
the boundary layer suggests a slightly different scaling 
relation for the eventual thickness h. If the boundary 
layer simply thickens diffusively prior to its collision with 
the top surface, then h should scale as V^^L^onset? where 
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vl is the kinematic viscosity of the hquid. Since tonset 
is simply a/U^ (see Fig. 2), we see that h/a should obey 
the Blasius scaling Re~^l^ . This is also consistent with 
our simulation results, basically because our range of Re 
is too limited to resolve a scaling exponent of 1/2 from 
the previously proposed 2/5 scaling (See supplementary 
materials). 

Before concluding, we comment on how these results 
on impact without air effects may relate to splash for- 
mation when air is present. In the simulated impacts, 
air effects are negligible and the boundary layer always 
remains attached to the wall. As t approaches 2r, the 
radial expansion slows and the pressure gradient within 
the liquid is essentially zero. In this time window, the 
boundary layer is not securely attached to the wall. Any 
external perturbation that adds an adverse pressure gra- 
dient, e.g. increased resistance from the air at larger am- 
bient pressures, may be enough to cause the outer edge 
of the boundary layer to separate from the wall. Since 
the surface profile is coupled to the boundary layer, the 
separating boundary layer may peel the thin liquid layer 
away from the wall, forming the beginning of a corona. 
Simulations to check this idea are underway [s^- 

In conclusion, we have simulated the impact of a vis- 
cous oil drop when the ambient air pressure is reduced to 
a very low value, so that impact at several m/s does not 
produce a splash. Results on the large-scale shape de- 
formation agree quantitatively with measurements from 
available experiments. The simulation reveals that the 
thin, spatially uniform pancake shape that a falling drop 
gets flattened into owes its existence to the boundary 
layer in the liquid drop created by impact. 
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